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SIGNAL GENERATION AND MIXING 
ELECTRONICS FOR FREQUENCY-DOMAIN 
LIFETIME AND SPECTRAL 
FLUOROMETRY 

The U.S. Government has a paid-up license to certain 
technologies disclosed in this invention and the right in 
limited circumstances to require the patent owner to license 
others on reasonable terms as provided for by the terms of 
the following contracts: NAS1-20426 and NAS1-20162 
awarded by NASA; DAAH01-91-R198 awarded by US 
Army Missile Command; DAAL06-92-C-0014 and 
DAAD07-91-C-0127 awarded by US Army White Sands 
Missile Range; DAAA15-93-C-0034 and DAAM01-94-C- 
0033 awarded by US Army Chemical and Biological 
Defense Agency; and F41624-95-C-6010 and F41624-97- 
C-6029 awarded by US Air Force Armstrong Laboratory. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention generally relates to fluorescence 
measurements. More specifically, the present invention 
relates to frequency domain measurements of the fluores- 
cence lifetime and the fluorescence spectrum. 

2. Description of the Related Art 

The process of fluorescence occurs when a substance, 
such as a molecule, absorbs light at one wavelength (or 
energy), and then emits light at a longer wavelength (or 
lower energy). A slight time delay occurs from when the 
substance absorbs light and when the substance re-emits 
light at the longer wavelength. This time delay is known as 
the fluorescence lifetime. 

FIG. 1 depicts the fluorescence process schematically for 
a diatomic molecule, which is a molecule that is comprised 
of two atoms. There are several "types" of energy that are 
internal to a diatomic molecule, one of which is vibrational 
energy. The classic example of vibrational energy is imag- 
ining the two atoms connected together by a spring where 
the atoms oscillate back and forth along the axis of the 
spring. As the two atoms approach each other, they experi- 
ence repulsion due to the proximity of their negatively 
charged electron clouds. As the atoms pull apart from one 
another, they experience an attraction that one can imagine 
as the result of the attractive forces between the positively 
charged nucleus of one atom for the negatively charged 
electron cloud of the other atom, and vice versa. The 
distance between the two atoms is known as the internuclear 
distance, which constantly changes as the two atoms oscil- 
late back and forth. As the internuclear distance decreases, 
a rise in potential energy occurs due to the repulsive forces; 
as the internuclear distance increases, a rise in potential 
energy again occurs due to the attractive forces. FIG. 1 
illustrates this difference in potential energy by line 11 in the 
potential energy "well" 10 and by line 13 in the potential 
energy "well" 12 where line 11 and 13 plot the potential 
energy as a function of internuclear distance. 

A group of horizontal lines 20, 22, and 24 appear in each 
of the potential energy wells 10 and 12. Each horizontal line 
represents the individual vibrational energy states possible 
for the diatomic molecule. Quantum mechanics requires that 
the vibrational frequency of the spring be within certain 
values. In other words, the spring may oscillate at frequency 
'a', and the spring may oscillate at frequency V, but it is a 
physical impossibility for the spring to oscillate at any 
frequency between 'a' and *b\ In FIG. 1, any two adjacent 
horizontal lines would represent 'a' and 'b\ 
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The different potential energy wells 10 and 12 represent 
another type of energy in the diatomic molecule: electronic 
energy. Different electronic energies occur when the mol- 
ecule absorbs energy in such a fashion that it causes an 

5 electron to move to a higher energy configuration within the 
molecule. The classical example of the energy absorption is 
the 'changing the spring' that connects the two atoms. If we 
"add" electronic energy to the potential energy well 10, the 
energy raises the potential energy well 10 to the energy level 

10 of potential energy well 12. 

With this background, we can use FIG. 1 to describe the 
process of fluorescence and fluorescence lifetime for a 
diatomic molecule. The example molecule is originally in a 
state depicted by the lower potential energy well 10. The 

1 5 molecule absorbs a photon of the correct amount of energy 
that induces an electron to move to a higher electronic 
energy state, represented by the upper potential energy well 
12. The diatomic molecule also tends to undergo a change 
from a ' lower vibrational state' (the horizontal lines 20 and 

20 22) in the lower potential well 10 to a higher vibrational state 
in the upper potential well 12. This initial excitation of the 
molecule is shown by line 14. Through any one of numerous 
possible processes, the higher vibrational state 24 in the 
upper well decays to a lower vibrational state 22 in the upper 

25 well (as shown by line 18). The process of non-radiative 
decay to a lower vibrational energy level in the upper well 
occurs very rapidly. After the non-radiative vibrational 
relaxation, the electron will want to return or revert back to 
a lower energy state as shown by line 16 (in other words, the 

30 electron wants to move back to where it was before the 
whole process started), and in so doing, the molecule emits 
a photon of light (with a lesser amount of energy than the 
excitation photon). This lowering of the energy state and the 
emission of a photon is the process of fluorescence. The 

35 process of emitting the photon has an associated time delay, 
which is the fluorescence lifetime. 

In practice, we typically probe many different types of 
molecules at once with the excitation light pulse. FIG. 2 
illustrates the case where a single short excitation pulse of 

40 light 30 is absorbed by a sample of identical molecules all 
at once. The fluorescence decay curve 32 resulting from a 
typical fluorescence response of a sample of identical mol- 
ecules is exponential in nature because not all of the iden- 
tical molecules emits its fluorescence photon at precisely the 

45 same time. The exponential decay follows a mathematical 
function so that we can calculate the fluorescence lifetime, 
x: 

50 /(/)=/ 0 -^- 5 ^ 

where ^ is the time the excitation pulse, I 0 is the initial 
fluorescence and I(t) is the observed fluorescence intensity 
as a function of time. 

55 In the cases with multiple molecular types (where more 
than one lifetime decay is present within a sample of a target 
mixture), the exponential decay seen in FIG. 2 would appear 
as a sum of exponential functions. Prior art fluorometers 
typically use a type of time-correlated (or time-resolved) 

60 system that count the emitted photons (from an excited 
molecule) in order to measure the fluorescence lifetime. 
Other prior art systems add the ability to take the fluores- 
cence lifetime measurements in the frequency domain by 
modulating the dynode of a photomultiplier tube, followed 

65 by a mixing and correlation procedure. These prior art 
systems are cumbersome, time consuming, and complicated 
to operate. The present invention overcomes the limitations 
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of the prior art systems by utilizing a novel technique to FIG. 1 illustrates the fluorescence process for a diatomic 

measure the fluorescence lifetime and spectrum. Instead of molecule. 

irradiating the target sample with a single short pulse of light pic 2 illustrates the absorption and emission of energy 

(photon counting), the present invention continuously irra- for the fl uoresce nce process. 

diates the target sample with a light source whose amplitude 5 mQ 3 iUustrat6S the chemometric extraction proceS s of 

modulation frequency is stepped with time. This technique c fl nn A i;f« 

allows us to u^e a chemometric analysis to automatically the P resent * ventl0n for florescence spectrums and life- 

extract the lifetimes from the 'phase delay ' and ' intensity vs. times. 

frequency' characteristics of the emitted light. FIG. 4 is a block diagram of the present invention. 

SUMMARY OF THE INVENTION 10 FIG - 5 k a block diagram of a component of the present 

m ... c , , invention for generating an excitation signal. 

The present invention is a system for chemometric analy- ° 0 

sis for the extraction of the individual fluorescence spectrum FIG. 6 is a block diagram of an embodiment of the present 

and fluorescence lifetime from a target mixture. The present invention for processing an emitted signal, 

invention comprises a processor with an apparatus for FIG. 7 is a block diagram of the present invention that 

generating an excitation signal to transmit at a target mixture illustrates the electrical pathways of the excitation signal and 

and an apparatus for detecting the emitted signal from the the emitted signal. 

target mixture. One embodiment of the present invention FIG. 8 is a block diagram of an embodiment of the present 

uses a processor that comprises a computer that extracts the invention for a chemometric analysis of the emitted signal to 

individual fluorescence spectrum and lifetime measurements determine the fluorescence spectrum and fluorescence life- 
from the frequency and wavelength data acquired from the 20 t ] mes 

emitted signal. The present invention first determines the G fig 9 is another block diagram of an embodiment of the 

and S matrices from the frequency and wavelength vectors. pfesent inveDtion for a chemometric analysis of the emitted 

A renormalization of w" S,. to G x occurs next due to the gi * tQ determine me fluorescence spectrum and fluores- 

diflerences in units. Next, the present invention determines c&Qce lifetimes 

the initial U, V, g and w^s approximations The present 25 ^ ^ { ability t0 dif . 

invention uses an iterative solution that first requires the ^ & ^ 

initialization of the decision variables including t). The r 

iterative solution compares the decision variables for con- DETAILED DESCRIPTION OF THE 
vergence to see if further approximation determinations of INVENTION 
the U, V, g and w _1 s are necessary. When the solution 30 ... ~ , . ■. 
conve ges, the present invention then determines the ne present invention is a system for chemometnc analy- 
reduced best fit error for the analysis of the individual » the extraction of the individual component fluores- 
fluorescence lifetime and the fluorescence spectrum. And «** *P«*a and fluorescence lifetimes from a target mix- 
finally, the system of the present invention extracts indi- tore. Additionally, the present invention is a method and 
vidual component fluorescence lifetimes and fluorescence 35 apparatus . for generating and mixing signals for frequency- 
spectra from the emitted signal of the target mixture. domain lifetime and spectral fluorometry. This disclosure 
r _ . ,JH. „ - , j +u a a describes numerous specific details that include specific 
The present mventon additionally includes a method and P ^ a 

apparatus for generating and mixing signab &r frequency- F > . o£ ^ For 

domain hfenmeandspectralfluorometry.ThepresenUnven- f S describes the extraction of 

SSSd Zt^^XSLtt the spectra and fluorescence hfeumes from a 

- f.i « . . , , K A target mixture. However, the practice of the present inven- 

frequency of the dnvmg/reference signal over time and a 7 s Tj .l ^ i • p 

. M . * , A . * *i. • • • ft1tkn f , tion includes the extraction or analysis of other physical 

signal generator mat generates the mixing ; ogoal that modu- * tfaan ^ 

lates the amplitude and/or the frequency of the mixing signal P J\ 0 ne skilled in the art will appreciate that one 
over time. The dereference ^ generator drrv«s an J^Jactice ^ invention ^ these ^ 
excitation signal generator that generates the excitation ~r/./ A ,,. + . ,f 4 .. , a^^a^^uZ^^ 
• rpi_ f- 4U % * ^ t™*™;** details. Additionally, this disclosure does not describe some 
signal. The present invention then directs and transmits the 39 T7 __. ^ k.*™*,— „. 
excitation signal at a target mixture, which absorbs the well known processe s such as Fourier ^foms 
energy from the excitation signal. The property of fluores- eigenequations, eigenfiinctions, eigenvalue^ ansfoms o 
cencf causes the target mixture to emit an emitted signal that * best * ***** m M m order not t0 obscure the P resent 
the present invention detects with a signal detector. The mV ^ n ° n ' . 4 
present invention uses a mixer that mixes the emitted signal ™ e P resent »7 nta » ^ 68 a ao ™ 1 s f tem mat ™£ a 
with the driving/reference signal to produce a processor chemometric : analysis to determine the fluorescence de- 
reference signal Another mixer mixes the emitted signal times and the fluorescence spectra of a target mixtai* 
with the mixing signal to produce a data signal. The present * Instead of irradiating the target sample with a single short 
invention then uses a processor to compL the processor of hght (photon counting) as other prior art systems, 
reference signal with the data signal by analyzing the the present invention continuously irradiates the target mix- 
differences in the phase and the differences in the amplitude ture with a hght source whose ^amplitude modulation fre- 
between the two signals. The processor then extracts the °W ■ ^pped ^h time. FIG. 3 literates the process of 
fluorescence lifetimes! fluorescence spectrum of the emit- « the present invention. The excitation hght or signal 42, 
ted signal from the phase and amplitude information using winch is * e ^ ^ P resent lnventlon ™* t0 c ? nt ^- 
a chemometric analysis ousl y irradiate the target mixture, is amplitude and phase 
' ' modulated in a sinusoidal fashion. As previously discussed, 
DESCRIPTION OF THE DRAWINGS the em i tte d signal 40 is the photon emitted (emission 
To further aid in understanding the invention, the attached 65 intensity) from the molecule as a result of the fluorescence 
drawings help illustrate specific features of the invention and process. The emitted signal 40 will appear as a delayed 
the following is a brief description of the attached drawings: sinusoidal-intensity light of a longer wavelength. Embedded 
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within the emitted signal is the information necessary to 
extract the fluorescence lifetime and spectrum. However, 
one must use a chemometric analysis to extract the fluores- 
cence lifetimes from the ' phase delay O' of the emitted light 
(compared to excitation signal) and the 'intensity vs. wave- 
length' (or spectrum) characteristics of the emitted light 
(compared to the excitation signal). 

FIG. 4 is a block diagram of the apparatus component of 
the present invention. A processor 50 couples to a data 
acquiring component 52 through a data input/output path 68. 
The preferred embodiment of the processor 50 comprises a 
computer. The data acquiring component comprises a gen- 
erator component 54 and a receiver component 56. Coupled 
to the generator component 54 is an excitation signal gen- 
erator 58. The excitation signal generator produces an exci- 
tation signal (or light or photon) that is transmitted to a target 
mixture 66. The target mixture 66 absorbs the electronic 
energy of the excitation signal as the molecules of the target 
mixture move to a higher electronic energy state. During the 
fluorescence process, the molecules of the target mixture 
emit photons (an emitted signal 64) as the molecules move 
to a lower electronic energy state. A detector 60 detects the 
emitted signal 64 from the target mixture 66. Coupled to the 
detector is the receiver component 56, which processes the 
returned emitted signal 64 into a reference signal and a data 
signal for the processor 50. The processor 50 then uses a 
chemometric analysis to determine the individual fluores- 
cence spectrum and fluorescence lifetime of one or more 
fluorophores within the target mixture 66. 

FIG. 5 is a block diagram of the generator component 54 
of the present invention for generating the excitation signal. 
The present invention uses heterodyning techniques to pro- 
duce two sinusoidal RF signals, a driving/reference signal 
90 and a mixing signal 92. The present invention modulates 
the frequency of the two signals from 10 MHz to 200 MHz. 
One skilled in the art will appreciate that the present 
invention could vary the signals over a much larger fre- 
quency range. The preferred embodiment of the present 
invention generates the two signals with a frequency differ- 
ence of 10 kHz. Another embodiment of the present inven- 
tion generates the two signals using an adjustable offset 
frequency where the offset frequency is set through to the 
present invention's control software. 

Referring back to FIG. 5, the present invention generates 
the driving/reference signal with a driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73. The driving/reference signal generator 
71 comprises an oscillator 74, a phase -locked-loop (PLL) 
70, a mixer 78, and a signal splitter 72. The mixing signal 
generator 73 comprises an oscillator 76, the phase-locked- 
loop (PLL) 70, a mixer 80, and the signal splitter 72. The 
present invention generates the two signals by mixing a 
signal from the PLL 70 with the signals from the two fixed 
oscillators 74 and 76. The preferred embodiment of the 
present invention operates oscillator 74 at 290 MHz and 
oscillator 76 at 290.01 MHz. One skilled in the art will 
appreciate that other operating frequencies for oscillators are 
possible. The preferred embodiment of the present invention 
can set the output of the PLL 70 to a frequency ranging from 
300 to 560 MHz, which the present invention controls 
through the control software of the present invention 
through the digital I/O datapath 68. 

Coupled to the driving/reference signal generator 71 is a 
filter 82 and an amp 86. The filter 82 filters out any unwanted 
components from the driving/reference signal, and the amp 
86 provides additional gain to the signal if necessary. 
Coupled to the mixing signal generator 73 is a filter 84 and 
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an amp 88. The filter 82 filters out any unwanted compo- 
nents from the mixing signal, and the amp 86 provides 
additional gain to the signal if necessary. The present 
invention uses the driving reference signal 90 as a phase and 
amplitude reference signal on the receiver component 56 (of 
FIG. 6). And, the present invention sends the mixing signal 
92 to the receiver component 56 where it is split. 

The present invention uses the driving/reference signal 90 
to directly drive and modulate an excitation signal generator 

) 58 (of FIG. 7) that may comprise a laser diode, a light 
emitting diode (LED), or a Deuterium lamp. The present 
invention can also use the modulation to directly modulate 
an electro-optic crystal for use in external modulation of gas 
laser sources or other continuous emitting excitation light 

; sources. 

FIG. 6 is block diagram of the receiver component 56 of 
the present invention for processing the returned emitted 
signal. The present invention receives the modulated emitted 
signal 64 from the target mixture through a detector (60 of 

) FIG. 7). The receiver component 56 receives the detected 
emitted signal 96 from the detector through an amplifier 97 
and receives the driving/reference signal 90 through an 
amplifier 99. The receiver component 56 additionally 
receives the mixing signal 92 through a splitter 98. The 

5 present invention then mixes the detected emitted signal 96 
with the mixing signal 92 with a mixer 102 and mixes the 
driving/reference signal 90 with the mixing signal 92 with a 
mixer 100. The resulting signals contain the sum and dif- 
ference frequencies of the mixer's (100 and 102) input 

3 signals. The present invention generates a data signal 114 at 
10 kHz by filtering the mixed detected emitted signal 
through a filter 106 with an amplifier 110 to eliminate the 
high frequency components. The present invention addition- 
ally generates a processor reference signal 112 at 10 kHz by 

5 filtering the mixed driving/reference signal through a filter 
104 with an amplifier 108 to eliminate the high frequency 
components. One skilled in the art after reading this disclo- 
sure will appreciate that one can use frequencies other than 
10 kHz to process signals. An analog to digital converter 118 

3 converts the data signal 114 into digital data for transmittal 
over the digital I/O path 68 to the processor (50 of FIG. 4), 
and an analog to digital converter 116 converts the processor 
reference signal 112 into digital data for transmittal over the 
digital I/O path 68 to the processor. The present invention 

5 then processes the two input signals to determine their 
relative amplitudes and phases because the amplitude and 
phase of the 10 kHz signals are directly proportional to the 
amplitude and phase of the high frequency modulation 
signals. As the modulation frequency is stepped from 10 

o MHz to 200 MHz, the processor reference signal 112 and the 
data signal 114 remain at 10 kHz. 

FIG. 7 is a block diagram of the present invention that 
illustrates the electrical pathways of the excitation signal and 
the emitted signal. The present invention generates the 

5 driving/reference signal with the driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73. The present invention uses the driving/ 
reference signal 90 to directly drive and modulate an exci- 
tation signal generator 58 that may comprise a laser diode, 

o a light emitting diode (LED), or a Deuterium lamp. The 
present invention can also use the modulation to directly 
modulate an electro-optic crystal for use in external modu- 
lation of gas laser sources or other continuous emitting 
excitation light sources. The excitation signal generator 58 

5 produces an excitation signal 62 that is transmitted to a 
target mixture 66. A photodetector 120 detects the excitation 
signal 62 on its way to the target and regenerates the 
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driving/reference signal 90 for retransmission to the receiv- line (linear least squares fit) to a graph of log(I(t)) vs time, 

ing component 56 of the present invention. The target we can determine the lifetime of a fluorophore. However, if 

mixture 66 absorbs the electronic energy of the excitation more than one fluorophore is present, the intensity variation 

signal as the molecules of the target mixture move to a with time depends on the concentrations, c„, and lifetimes, 

higher electronic energy state. During the fluorescence 5 V of tne individual fluorophores: 
process, the molecules of the target mixture emit the emitted 

signal 64 as the molecules move to a lower electronic energy /(/> = c„«- (, -^ ,/T « ( 4 ) 

state. A detector 60 detects the emitted signal 64 from the * 
target mixture 66. The detector 60 could comprise a photo- 

detector. Additionally wavelength dispersion of the emitted 10 Nqw ^ d ^ must be & b a QQn ^ m3LT 

signal can be accomplished by using a spectrophotometer, least es ss? which requires ioitial esses for the 

monochromator or low-pass, high-pass and notch optical number of fluorophores? r , the ir concentrations, and their 

filters. From the detector 60, the present invention mixes the lifetimes. In general, convergence of the solution cannot be 

signals as previously described for FIG. 6 above. The guaranteed as it depends on the accuracy of the initial 

present invention then uses a chemometric analysis to deter- 15 guesses. 

mine the individual fluorescence spectrum and fluorescence Weber showed that decay constants and relative intensi- 

lifetime of one or more fluorophores within the target ties of an arbitrary number of independent components in a 

mixture 66. heterogeneous fluorescence emission could, in principle, be 

The present invention comprises several techniques or determined by measuring phase shift and relative modula- 

alternative embodiments to "target" the target mixture. One 20 tion of the total fluorescence at an equal number of frequen- 

embodiment of the present invention uses a target compart- cies. See Weber, G., Resolution of the Fluorescence Life- 

ment for holding the target mixture within a chamber. A times in a Heterogeneous System by Phase and Modulation 

container may contain the target mixture, or the mixture may Measurements, J, Phys, Ckem., 85, 949-953 (1981). At each 

be free floating, in a fluid state, within the chamber. Another frequency, the phase and modulation measurements yield 

embodiment of the present invention uses an optical probe 25 the complex Fourier transform of the fluorescence impulse 

that combines the signal generator 58 and the detector 60 response, G+iS. The moments of a distribution of lifetimes 

into a single probe. A user of the present invention would are found as linear combinations of the G's and S's and 

then probe or place the optical probe into the target mixture. Prony's method is used to obtain lifetimes and concentra- 

And, another embodiment of the present invention may use tions of the components from the moments. Due to its 

a laser to target a free floating target mixture, such as the 30 sensitivity to noise, this method is really not practicable for 

emission from a car (at a distance). One skilled in the art will more than binary mixtures. 

appreciate that other embodiments of targeting the target Lackowicz and his coworkers developed an analysis of 

mixture are possible after reviewing the disclosure of the frequency domain phase and modulation data based on a 

present invention. weighted non-linear least square analysis. See Lackowicz, J. 

FIG. 8 is a block diagram of an embodiment of the present 35 R > et aL > Analysis of Fluorescence Decay Kinetics from 

invention for a chemometric analysis of the emitted signal to Variable-Frequency Phase Shift and Modulation Data, Bio- 

determine the fluorescence spectrum and fluorescence life- P h y s - J > 46 > 463-477 (1984). By taking data at several 

times. Before discussing the block diagram of FIG. 8, we emission wavelengths, they were able to resolve binary 

first need to discuss the mathematical foundation used for mixtures whose components differed in lifetime by 30% and 

the chemometric analysis. Fluorescence occurs when a 40 teraar y mixtures the ratio of whose range of lifetimes 

molecule, excited by a photon to a higher electronic state, spanned a decade. Best fit is determined by minimizing a 

decays non-radiatively to a lower vibrational level of the reduced chi-squared function. Poor initial guesses can still 

excited state, and then decays radiatively to a vibrational produce unpredictable results. 

level of its ground electronic state. The non-radiative decay Working with an Excitation-Emission-Frequency Array, 
is assumed to be 'fast'. The characterization of the radiative 45 McGown and her coworkers employed a principal compo- 
decay is the fluorescence lifetime, x. And as previously «ents analysis followed by a non-linear least squares mini- 
stated, the wavelength of the fluorescent photon is greater mization on an over determined set of data to extract the 
(meaning less energy) than that of the exciting photon. s P ectra and lifetimes of components in a mixture. See 
If we apply a 'monochromatic', short light pulse to the c Burdi * ] etal Resolution of Mumcomponent Fluo- 

molecules of a fluorescent material at a time to, ideally 50 ™ cent M ™ b J ^ S1S of the ^f^™" 
represented by Frequency Array, J. Chemometncs, 4, 15-28,(1990). In 

* addition to requiring a non-linear minimization, a disadvan- 

tage of this method is that it requires evaluation over a large 
E{t)=E 0 S(f-f 0 ) (l) (3 dimensional) parameter space, 

where E 0 is the total energy in the pulse and E(t) is the 55 ^ P resent invention utilizes a novel system that uses a 
intensity of the incident pulse as a function of time, we find chemometric analysis for measuring the fluorescence life- 
that the fluorescent intensity at a fixed wavelength has the time and tne fluorescence spectrum of a target mixture. 
£ orm: Instead of irradiating the target sample with a single short 

pulse of light (photon counting) as other prior art systems, 
60 the present invention continuously irradiates the target mix- 
/(rHo-e" 0 "' 0 ^ ( 2 ) ture with a light source whose amplitude (intensity) and 

where phase (frequency) are modulated (varied) with time. In 

evaluating experimental phase resolved fluorescence spec- 
troscopy (PRFS) data, the variables sought are the 
T <* =cvE ° ( 3 ) 65 concentrations, c„ the lifetimes, x t , and the emission spectra, 

and a has appropriate units. Thus, by measuring the decay ot/X), of the individual fluorescent contributors. The vari- 
of the fluorescent intensity with time and fitting a straight ables that we may control include the emission wavelength, 
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X, the modulation frequency, oo, and the detector reference From which we see that G and S are proportional to the 

phase angle We first consider an amplitude modulated Fourier transform of Equation (11). 

excitation: In the case of noiseless data, where r fluorophores are 

present, it is necessary to make measurements at r values of 

5 X and r values of w to obtain: 

E(£)=EJ[1+M e sin (<at)l t>0 (5) 

where oa is the modulation angular frequency, M e is the g=cc-a and s=cs-A, (12) 
modulation amplitude and E 0 is the average amplitude of the 

excitation. where 

At steady state, the modulated part of the fluorescent 10 
intensity may be given by 

/xO^oM^A sin (<0f-$) (6) HS jk \ 

where A is the amplitude response at the fluorescent wave- 15 cc=[cos 2 ^^HccJ 

length and <D is the phase lag at this wavelength and cj=[cos sln ^)H«J 
frequency. 

We originally detected this signal (the emitted signal) by A={e t (X k )Ha tk \ (13) 

correlating it with a reference signal (the driving/reference . . . 

signal) bfmixing the two signals with the mixing signal and 20 where the quantities defined in boldface are matrices. The 

integrating them over some multiple of a period. The cor- matrix indices i, j, and k run over the components, the 

related signal could therefore be expressed as: modulation frequencies, and the emission wavelengths, 

respectively. 

2k If no noise is present, we may take the indices i, j, and k 

S\{a>, $r) = — P 7 /*(r)sinM -* R )dt 25 above to run from 1 to r, the number of components. As the 

27 Jo fluorophores are assumed to be independent and as the 

= E 0 M e A<x*t$ R - $) frequencies {coy} may be chosen arbitrarily, all the matrices 

defined above have true inverses. Thus, from Eqns. (12) and 

= \e 0 M c y&cos^costo - (13) we find: 
2 it 30 



cc^-G^A (14) 



S=cs-cc- x 'G 



where r is the number of components present (or 
fluorophores) and ® R is the phase of the reference signal and 
relative to the backscattered excitation signal. We observe 35 
that S x (co,&^) may be written as: 

1 (8) 

40 

S'G~ 1 'cc=cs 

Therefore, as an alternative to Eqn. (8), we may determine 
G(co,X) and S(co,X) to within a common normalization factor, n 0 w, we note that: 
E 0 MJ2, for any fixed w and X by making measurements at 

two different, but known, values of ® R . In the presence of [cos^w^sin^ (15) 

noise, more reference phase values would be used in the 45 

measurements and fit in a least squares manner to G and S. = [cosV t (<^)tan^(w,) 

We then find that for each value of co and X (emission 
wavelength) we have (modulo the normalization factor): gy from m we have; 

ChD^^M W 50 cs=[(coapM ^ l0j r l ] (16) 

5(M,A) = 2e,a)cos^,( M )sii^,(«) <10) 

55 where w»diag[co ; -] and T«diag[t-]. Hence, we find 

with tan ^~(xrt £y where i runs over the number of fluoro- 
phores present and assumed to be acting in an independent 
(uncorrelated) fashion. 

The preferred embodiment of the present invention allows ^ or 
us to determine G and S by mixing the emitted signal from 
the mixture down to an intermediate frequency of 10 kHz, 

sample and digitize it at 50 kHz, Fourier transforming the w-^S-G^-cc^cc-T 

result and picking off the signal at 10 kHz. When this is ^ q . g . q ^ fonn of a standard eigenvalue/ 

done the mixed signal has the form: &5 eigenfimction equatkm5 where ^ X( m the eigenva i ues 

the columns of cc are the eigenfunctions. We may also write 
f x (t)=EMG sui (<* JF tys cos («w)] (11) this in the form: 



^==[<cos^(w,)> 2 <cJ,rJ 
= w-cc-T 
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number of components. Furthermore, some degree of noise 
will be present in the experimental data. Therefore, it is 
cc-^w-^s-cr^cc^ (18) useful to determine a means of solving for lifetimes and 

spectra when the number of frequencies and the number of 
emission wavelengths are unequal and when each is greater 
than the number of independent components present in the 
mixture. 

In the following, we shall therefore assume that r com- 
^ ponents are present in a mixture, for which we have taken 

10 data at N frequencies and M emission wavelengths. We 
(20) assume that N,M^r and that N*M, in general. 
Tr(w~ l -S- cr l ) = £ T| From Eqns. (12) and (16), we therefore have: 



From which we have: 

r 



(19) 



Since c^/coy) is known once x t - is known, we may determine 15 



G=ccA (28) 



cc or cs and solve the equation: (NxM)-(Nxr) (rxM) 

and 

G=cc-A or S=cs-A (21) 

20 

for A with a standard LU decomposition, which gives us s=wcc-TA, (29) 



(NxM)-(NxN)'(Nxi)-(rxr)-(rxM) 



A^-ccprw^S (22) 



where the number of rows and columns of each matrix is 
where ccpi is the pseudo inverse of cc. 25 indicated below it. In the noise free case, all the matrices are 

Another, equivalent, way of writing the eigenfunction of rank r. In the presence of noise, G and S will be of rank 
equations proceeds from the observation that: N or M, whichever is smaller. 

In the absence of noise, all matrices are of rank r. We 
manipulate Eqns. (28) and (29) slightly, to write: 

GA- X =cc (23) 30 

and G-ccIA 



and 



S=csA=wcc-T-A (24) 



35 



so that W^S-cc TA (30) 

where I is the rxr identity matrix. 
w-i-S-A- 1 =cc-T=<j'A- t 'T (25) We may now make a singular value decomposition on the 

40 first of Eqns. (30): 



or 



G-^w^S-A-^-^T (26) 



G=U-C 1 -V T (31) 
(NxM)-(Nxr) (rxr)-(rxM) 



Solving this eigenequation determines x t as the eigenval- 45 
ues and A" 1 as the eigenvectors. However, as the eigenvec- where c i * the ^agonal matrix of singular values of G and 
tors are normalized, A may be better determined as above if u and v are eacn column orthonormal, i.e.: 
information about the component concentrations is required. 

Once the spectrum and lifetime of each component is u T -U-v r 'V-l (32) 

known, we may use the information to identify the fluoro- 50 

phore. The concentration of each fluorophore may then be From Eqns. (31) and (32), we have: 
determined by comparing its spectrum against its standard 
spectrum. We find: 



^ _ g t (A)/<MA)r, (27) 
Cj ej(X)/aj{\)Tj 

with 



55 



C^U^G-V (33) 

We define: 

c^iF-w-^s-v (34) 



60 where C 2 is not diagonal, in general, but is rxr. 

It may easily be shown, by a process that is similar to 
Gram-Schmidt orthogonalization, that while, in general, 
where a t {\) is the standard spectral value of the ith com- neither matrix ^ invertiblej cc has a left i^se and A has 
ponent at wavelength X and x i is its lifetime. a inverse* 

The number of components in a mixture may not be 65 
known a priori, and the method requires the number of 

frequencies and wavelengths employed to at least equal the cdi-cc=i and AAn=i (35) 



cc^U-Py and A=P V -V T (36) 
with 5 
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From Eqns. (30) and (31), we may therefore write: Asecond method for solving Eqns. (30) is to use Eqn. (31) 

to write the pseudo inverse of G: 

Gpi^V'C^-U 7 (43) 

so that we have: 

P^C^Ari and Py=cclrU-C v 

where Y v and P v are each rxr. G-Gj^and GpiG-vv', 

We may interpret the first of Eqns. (36) in terms of its 10 where U-U r is and NxN matrix and V- V r is an MxM matrix, 
column vectors as the representation of the r vectors of cc in an d eac h c f t h ese matrices is of rank r. 
the orthonormal basis of the r vectors of U: From Eqns. (30) and (35), we may write: 

^ (37) « 

£,=2j^A 15 GAri=cc (44) 

Substituting this into the first of Eqns. (30), we have: 

where c y is the jth column of cc, U, is the ith column of U 

and P t y is the corresponding element of P^. From this Gpt-w-^s-An^v-v^AH-T (45) 

expression and the linear independence of the c, and U £ , we 

see that P,, is invertible. Linear independence of the c, From ^ ( 30 > and ( 36 ) we have: 
requires that: 

iv' 1 S^U-Pu-T-Py-T. (46) 

(38) 25 

2^=0 ™ l y if a j = Q all /=l,...,r. From Eqn. (32), we therefore find: 



From Eqn. (36), we see that Eqn. (37) implies: 



w - 1 'S=w- 1 -S-VV r . 

30 Substituting this into Eqn. (45), we therefore have: 



Gpi-w'-S'V-V T -Ari=V-V T -AH-T (47) 

Eqn. (47) is the eigenvalue equation for the matrix, 
From the linear independence of the U f , we must have 35 Gpi-w^-S, with eigenvalues equal to the diagonal elements 

of T and eigenvectors proportional to the columns of 
' (39) V-V r *Ari. The spectra, A, may be found as described above. 

Zj p ij a j = °» for all i - 1, r. When determining the lifetimes, we can renormalize Eqn. 

(47) by using 



j=i 



Thus, if the condition of Eqn (37) is to be satisfied, the pseudoinvG?) w^s evec=evec T (48) 
system of Eqns. (39) must have only the trivial solution. This 

is equivalent to the requirement that where evec=VV r Ari for T=diag(x 1? . . . , x r ) 

For the case when noise is present, let us assume that r 

45 components are actually present, and let us take N<M. The 

det [ p t/^°- experimentally determined matrices, Gx and Sx, are there- 

Therefore P^ has an inverse. A similar argument shows fore of rank N - We ^ t0 determine the NxM matrices, G 

that P is also invertible anc * ^» °^ ran ^ r rmrnrn i zc some suitably defined error. 

From Eqns. (30), (33)' and (34), we have: <n . We shaU define the Euclidean norm of a matrix, M, to be 

50 given by: 



C^P^T-Py and C^P^PPy. (40) 

Therefore, we may write: 

C^C^Py^TPy (41) 

or 

C-^Py^Py^T. (42) 

Eqn. (42) is the eigenvalue equation for the matrix, 
C^-C,, with eigenvalues equal to the diagonal elements of 
T, and eigenvectors proportional to the column vectors of 
P F \ The spectrum, A, may be determined from Eqns. (13) 
and (30) since cc is determined once the lifetimes are 
known. 



||M|i 2 = Tr{M T ■ M) = 7>(M - M T ) = £ £ |m„| 2 (49) 
» j 

55 

The error we define will depend on our knowledge of the 
statistics of the experimental data. We shall assume that the 
errors at each data point are uncorrelated with those at other 
data points and that the standard deviations from the sample 

60 means are the same for all data points. This is tantamount to 
the assumption that the statistics are independent of the 
modulation frequency and the emission wavelength for the 
ranges of the variables explored. Previous experience with 
experimental data tends to bear this assumption out. 

65 Given these assumptions, it is reasonable to define the 
error, / 2 , in terms of the equally weighted distance of the 
experimental points from their analytical counterparts, i.e.: 
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and M^x^QuGx+iw-^SxY'Qv-iw^-Sx) 



(66) 



(50) 



where G satisfies Eqn. (31) and where, from Eqns. (40) and 
(46), we see that S satisfies: 



w~ 1 -S=U-Ci'V r 



(51) 



We note that M v is an NxN matrix, M^is an MxM matrix, 
and both matrices are symmetric. 

If we assume that V is known, we may determine U by 
requiring that it be chosen to maximize Eqn. (63) subject to 
the constraint imposed by Eqn. (32). We may solve for U by 
means of Lagrange multipliers: 



Define: r\'*r)-T?[L v iU T -U~r)] 



(67) 



xHGx-u-c^vYMl^-sx-u-CzV 1 )? 



Therefore, from Eqns. (31), (50) and (51), we must find U, ™ 
V, C ± and Q, to minimize: 

where is a symmetric rxr matrix formed from the 
multipliers. 

(52) We may now extremize if with U unconstrained to find: 

15 

Minimizing Eqn. (52) will yield a least squares fit for the 
data. M v -u=U'L v (68) 

For Eqn. (52) to be a minimum, its partial derivatives with 
respect to each free variable must vanish. Applying this Since ix » symmetric, may be diagonalized by a 
requirement with respect to the variables of C 19 we find: 20 simil anty transform with an appropriate orthogonal matrix. 



or, from Eqn. (31), 

G=U-U T 'GxW T =Q u -G*Q v 

where 

Q ir U-U T ^d Qy=VV r 

Proceeding similarly for the variables of C 2 , we find: 

or w-^S^Qy-w-^Sx'Qv 
or ^U-lf-w-^S'Vx V 7 



(53) 



Therefore, we may write: 



(54) 



(69) 

25 where D is an rxr diagonal matrix and O is an rxr orthogonal 
matrix. 
We may then write: 



30 



M v -UO=U-0'D 



(70) 



(55) 
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(56) 
(57) 

(58) 40 



Substituting Eqn. (70) into Eqn. (63) we find: 

7j = Tr(MvU-0-0- 1 -U T ) 

= Tr[(U'0)'D.(U-0) T ] 
= Tr[D-{U-0)-(U-0)] 



or, finally, 



Eqns. (54) and (57) determine C 2 and once U and V are 
determined. We now proceed to determine these variables. 

By substituting for G and w -1 S from Eqns. (54) and (57) 
into Eqn. (50), we find that minimizing that expression 
implies that we must maximize: 



r\-Tr(D). 



(70) 



45 



r]-\\Q u -Gx'Qvf + \\Qu^- 1 Sx-Qyf 

or 

Observing that: 

Qv=Qu T , QuQu=Qu> 

and Qv=Qv T , Qv'Qv=Qv> 

we may rewrite Eqn. (59) in the alternative forms: 

r)=Tr(My<U-U T ) 
or v^TriMu-V-V 1 ) 

where M v ^x-QyGx r +(w- i -SxyQy(w- 1 -Sx) T 



(59) 



Thus, to maximize tj we choose U to maximize D as 
determined by Eqn. (70). We may easily do this by observing 
that Eqn. (70) has the form of an eigenvalue equation. The 
matrix D yields r of the N possible eignevalues of the matrix 
M^. We therefore solve the eigenvalue equation for M v : 



50 



MyU=U-E 



(71) 



(60) 



(61) 
(62) 



(63) 
(64) 
(65) 



taking the largest r eigenvalues of E for D and, without loss 
of generality, choosing U to equal the corresponding eigen- 
vectors of U. 

55 A similar procedure may be used to determine V from 
Eqn. (64) if U is known using the relationship M^ V^V-L^ 
to find L^. 

FIG. 8 is a block diagram of the preferred embodiment of 
the present invention that uses the above foundation as a 

60 basis for a chemometric analysis of the emitted signal to 
determine the fluorescence spectrum and fluorescence life- 
times. The present invention first acquires the data from the 
data signal and the reference signal to determine the fre- 
quency and wavelength vectors 132, Next, the present 

65 invention assumes an arbitrary but fixed number of 
fluorophores, r, 134. We then determine the initial G and S 
matrices 136. 



5,981,957 

17 18 

As the frequencies in w are large, we renormalize Since we collect an amplitude and a phase value for each 

||w- 1 .Sx|| to equal ||Gx|| so as to avoid skewing the results in emission wavelength and each modulation frequency, we 

favor of the Gx data 138. This amounts to a resealing of the have: 
units for the frequencies and the lifetimes. This resealing of 

units is compensated once the lifetimes are found, so they 5 

are expressed in seconds. tf*"2/VV x ( 76 > 

We begin the iterative approach by determining the initial where N^is the number of frequencies and N x is the number 

approximations for U, V, g and w~ 1 s approximations 140, of wavelengths used. 

142, and 144. The initial iteration begins by choosing the r From Eqns. (54), (57) and (68), we see that the free 

columns of U 1 to be the eigenvectors corresponding to the 10 parameters being fit to the experimental data are the ele- 

r largest eigenvalues of the matrix men ts of U, V, Ly and L u subject to the constraints of Eqns. 

(32). For r components, we therefore find that: 

Gx-Gx^w-^SxXw^-Sx) 7 " (73) 

or [G-G^w^-SHw^-S^-^UL (74) l5 N f = r ' N w +r-N k +2^-r(r+l)-2.^.r(r+l) (??) 

where L=diag(X, a , . . . , X r ) for \j> . . . >k r . =r-N w +r-N x 
Next, U 1 is used to form the matrix M m in accordance with 
Eqn. (66). We can determine V 2 from M m by solving the 

eigenfunction equation M^-V^-L^ by choosing its r col- We evaluate y£ R for each r, and choose the value of r 
umns to be the eigenvectors corresponding to the r largest 20 corresponding to the smallest value of x 2 * to represent the 
eigenvalues of M^. We then form M w in accordance with number of components which best fits the data 158. 
Eqn. (65). If two or more lifetimes are the same, their components 

We next initialize the decision variables 146 to produce will not be separable by this method, nor will they be 
the solution in a reasonable time, with a reasonable accuracy, separable by the method of PREEMs. Their contributions 
and a reasonable number of iterations. We begin the iterative 25 will remain combined as that of a single component with the 
solution by comparing the decision variables 148 for con- common lifetime and a spectrum given by the combined 
vergence to see if further determinations for U, V, g and w _1 s spectrum appropriate to their concentrations, 
are necessary as above. In the iterative solution, we proceed FIG. 9 illustrates an embodiment of the present invention 
to maximize t\, given by Eqn. (59), and recognizing that: for the method described in FIG. 8 for a chemometric 

30 analysis of the emitted signal to determine the fluorescence 
o^HGxlp-Hjw^Sxil^N (72) spectra md fluorescence lifetimes. The user interface 170 

allows the user of the present invention to interact with the 

Proceeding in this way, we form a sequence of rj values system of the present invention. The user interface allows 
{%} • We see easily that the sequence is monotonic non- the user to choose the number of components or fluoro- 
decreasing. By Eqn. (72) and the Bolzano-Weierstrass 35 phores to solve for. The data acquisition component 172 uses 
theorem, the sequence has at least one limit point in the the processor 50 coupled to a data gathering component 
interval [0,N]. The monotonic nature of the sequence guar- through a data input/output path 68 to acquire the data from 
antees that there can be but one limit point. Therefore the the emitted signal to determine the frequency and wave- 
sequence must converge. length vectors. A determiner 174 determines the G and S 

The values of U and V corresponding to the limit of the 40 matrices from the frequency and wavelength vectors. A 
sequence are then used to form G and w _1 S in accordance renormalizer 176 then renormalizes w"^ to G^. A deter- 
with Eqns. (54) and (57). Eqn. (47) may subsequently be miner 178 next determines the initial U approximation, and 
solved for the lifetimes and spectra as described previously a determiner 180 determines the initial V approximation, 
160, 162, and 164. and a determiner 182 determines the initial g and w" 3 s 

In the preceding, we assumed an arbitrary but fixed value, 45 approximations. An initializer/updater 183 then initializes 
r, for the number of fluorescing components present. Fits to the decision variables. The iterative solution of the present 
the data may be found for different values of r. The question invention uses a comparator 184 that compares the decision 
then arises as to the value of r that best fits the experimental variables for convergence to see if further determinations for 
data. Two qualitative criteria may be stated: (1) the spectra U, V, g and w^s are necessary. If the solution does not 
should look reasonable; and (2) the lifetimes should be 50 converge, the present invention uses the determiner 178 to 
positive real numbers. determine the next U approximation, the determiner 180 to 

Somewhat more quantitatively, we recognize that % 2 , determine the next V approximation, and the determiner 182 
suitably renormalized for frequency as described above, to determine the next g and w _1 s approximations. The 
represents our error for a fixed value of r. However, the value present invention then uses the initializer/updater 183 to 
of x 2 depends on the number of free parameters available for 55 update the decision variables. 

the fit, which, in turn, depends on r. Instead, we attempt to When the solution converges, a determiner 186 deter- 
remove this dependence by dividing by the difference mines the reduced best fit error for the analysis of the 
between the number of experimental data points and the individual fluorescence lifetimes and fluorescence spectra, 
number of free parameters present: The present invention then uses an extractor 188 that 

60 extracts the individual fluorescence lifetimes from the target 
y2 _ (75) mixture, and an extractor 190 that extracts the individual 

N x - N f fluorescence spectra. 

FIG. 10 illustrates the present invention's ability to iden- 
tify and discriminate between individual overlapping spec- 
where % R is the reduced error term, Ny is the number of 65 tral components in a target mixture. The present invention 
experimental data points and N f is the number of free allows for the extraction of both the individual fluorescence 
parameters. spectra and lifetimes from a target mixture of fluorophores. 



